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'Island Surfing' Mechanism of Electron Acceleration 
During Magnetic Reconnection 

M. Oka, 1 M. Fujimoto, 2 I. Shinohara, 2 T. D. Phan 1 

Abstract. One of the key unresolved problems in the study of space plasmas is to ex- 
plain the production of energetic electrons as magnetic field lines 'reconnect' and release 
energy in a exposive manner. Recent observations suggest possible roles played by small 
scale magnetic islands in the reconnection region, but their precise roles and the exact 
mechanism of electron energization have remained unclear. Here we show that secondary 
islands generated in the reconnection region are indeed efficient electron accelerators. We 
found that, when electrons are trapped inside the islands, they are energized continu- 
ously by the reconnection electric field prevalent in the reconnection diffusion region. The 
size and the propagation speed of the secondary islands are similar to those of islands 
observed in the magnetotail containing energertic electrons. 



1. Introduction 

The origin of energetic electrons of up to MeV in the 
Earth's magnetotail has been an outstanding problem for 
decades [e.g. Sarris et al., 1976; Terasawa and Nishida, 1976]. 
Although many theories consider magnetic reconnection as 
the mechanism of explosive energy release phenomena in the 
Earth's magnet osphere, the exact mechanism by which high 
energy electrons are produced is still unclear. 

A statistical survey by the Geotail spacecraft shows that 
accelerated and/or heated electrons are more likely to be 
detected around the outflow regions somewhat away from 
the center of magnetic reconnection [Imada et al., 2005]. It 
was suggested that electrons pre-energized at the X-lines 
further gain energy in the flux pile-up region especially in 
the earthward side of the X-line. While a Cluster event at 
the near-Earth reconnection site also confirms this scenario 
[Imada et al., 2007], another event by the Wind spacecraft in 
the distant magnetotail observed increasing flux of energetic 
electrons with decreasing distance from the X-line, with no 
clear evidence for a flux pile- up region in the outflow jet, sug- 
gesting that the region around the X-line is the dominant 
source of energetic electrons [0ieroset et al., 2002]. The dis- 
crepancy between the two events is not well understood so 
far. 

On the other hand, recent measurements by the Clus- 
ter spacecraft pointed out that magnetic islands may be an 
important agent responsible for the production of energetic 
electrons [Chen et al., 2008; Retino et al., 2008]. The main 
features of the observed islands are that they are small - of 
the order of ion inertia length - and that they move fast, 
close to the Alfven speed. These features suggest that they 
are the magnetic islands that are naturally formed as a con- 
sequence of secondary tearing instability of a thin current 
sheet [Eastwood et al., 2007, and references therein]. 

As for theoretical models of electron acceleration within 
magnetic islands, a contracting motion of volume filling 
magnetic islands has been considered [Drake et al., 2006]. 
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Because of the contracting motion, inductive electric fields 
are created at each end of the islands and particles that are 
fast enough to circulate inside the islands can receive repeti- 
tive 'kicks' from the electric fields. However, the presence of 
volume filling islands and their contracting motion are diffi- 
cult to verify observationally and have not been reported. 

In this Paper, we present an alternative scenario of elec- 
tron acceleration by performing 2D particle-in-cell (PIC) 
simulations of magnetic reconnection. We found that elec- 
trons are trapped in a secondary magnetic island so that 
they can be continuously energized by the reconnection elec- 
tric field in the diffusion region. 

2. Simulation 

We utilize a two and half dimensional, fully relativistic 
PIC code [Hoshino, 1987; Shinohara et al., 2001]. The ini- 
tial condition consists of two Harris current sheets. The 
anti-parallel magnetic field is given by B y /Bo = tanh((x — 
xr)/D) — tanh((x — xl)/D) — 1 where Bo is the magnetic 
field at the inflow boundary, D is the half-thickness of the 
current sheet and L x and L y are the domain size in x 
and y direction, respectively. /4 and x R =3L x /4 are 

the x-positions of the left and right current sheet, respec- 
tively. Periodic boundaries are used in both directions. The 
ion inertial length di is resolved by 25 cells. D=0.5di and 
L x = L y — 102 Adi. The inflow, background plasma has the 
uniform density of Nbo=0.2No where No is the density at 
the current sheet center. The ion to electron temperature 
ratio is set to be Ti/T e —b for the current sheet and Ti/T e =l 
for the background. The background to current sheet tem- 
perature ratio Tbk/T cs =0.1. The frequency ratio cj pe /^ce=3 
where oj pe and Q C e are the electron plasma frequency and the 
electron cyclotron frequency, respectively. The ion to elec- 
tron mass ratio mi/m e —2b and the light speed c is 15V a 
where Va is the Alfven speed denned as Bo / ^/A.nNomi. We 
used average of 64 particles in each cell. 297 particles per 
cell represents the unit density. 

No magnetic field perturbation is imposed at the begin- 
ning so that the system evolves from a tearing mode instabil- 
ity due to particle noise. Note that conventional simulations 
initiate magnetic reconnection by a small magnetic field per- 
turbation in order to save computational time [e.g. Birn et 
al., 2001]. Our simulation setup does not use such trigger so 
that the system does not produce any large initial velocity 
flows resulting from a lack of pressure balance. 

Figure la shows a snapshot of the out-of-plane component 
of electron current density J e , z obtained when the system is 
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Figure 1. (a) The out-of-plane component of the 
electron current density at Q c it=115. Superposed on 
the image are the locations of the most energetic elec- 
trons (e >1.4m e c 2 ). (b) Electron energy spectra in- 
tegrated over the rectangular regions indicated in (a) 
(solid curves). The dashed curve shows the spectrum at 
a square region centered at (x, ?/) = (25.6di, 26.4di) with 
each side 1.6di. 



well developed (Q c it=115). The diffusion regions are rep- 
resented by the localized current densities in a thin current 
sheet, and embedded within the diffusion regions are the 
small scale secondary islands (highlighted by the arrows). 
Superposed on the plot are the locations of the most ener- 
getic electrons (e >1.4m e c 2 , indicated by the white circles). 
Most of the energetic electrons are located within the large 
islands but some electrons do exist in the small, secondary 
island of the left current sheet. 

Those particles within the large islands are energized 
through the so-called multi-island coalescence resulting from 
the non-linear evolution of the tearing mode instability. A 
variety of different energization mechanisms has been found 
and they are fully described in our companion paper [Oka 
et al., 2010]. On the other hand, the energization of par- 
ticles within the small islands are not the consequence of 
multi-island coalescence and the physics involved are en- 
tirely different, as will be shown below. Hereafter, we focus 
on the electron energization process within the small scale, 
secondary islands. 

Figure lb shows the electron energy spectra obtained at 
the time of Figure la. It is evident that, at Q ci t=l 15, the 
highest energy electrons in the diffusion region are all con- 
fined in the small magnetic island, particularly at the very 
center of the island (red curves). This island overlaps the 
electron-scale diffusion region (or the so-called 'inner diffu- 
sion region') where electron acceleration should be efficient. 



For a reference, we also plotted a spectrum for the immedi- 
ate upstream (blue curve) and the downstream of this region 
(green curve). 

In order to clarify the history of electron acceleration, we 
followed the trajectory of the highest energy electron in the 
rectangular box in Figure 1. In Figure 2, we divided the tra- 
jectory into five segments and plotted over the image that 
shows the out-of-plane component of the electron current 
density. Note that the simulation itself is two-dimensional 
in x and y but the displacement Az can be calculated by 
integrating v z over time. The background images are the 
snapshots taken from each segment so care should be taken 
when comparing the trajectory with the images. The y po- 
sition as a function of energy e is also shown in the right 
hand panel. 

The secondary magnetic island starts to appear at 
Q c it^89 (Figure 2a). The electron acceleration seen up until 
Q c it~89 is what is typically seen in the inner diffusion re- 
gion. Had the electron left the diffusion region upward soon 
after Q c it^89, it would not have become one of the most 
energetic electrons in the system. In reality, it was deflected 
downward back to the diffusion region by the Lorentz force 
due to magnetic field of the emerging island. Then, a clear 
indication of the interaction between the electron and the 
secondary island is seen at Q c it^92 (Figure 2b). A distinct 
change of the electron behavior starts at Q c it~95. The elec- 
tron is trapped in a well developed island (Figure 2c). The 
island continues to develop so that the magnetic island be- 
comes as large as ^5di although electron current is rather 
localized of the order of 2di (Figure 2d). During this later 
phase, the electron continues to gain energy constantly (Fig- 
ure 2f). 

The behavior of the electron can also be studied from the 
time profiles of the electromagnetic field felt by the particle 
(Figure 3b). Before being trapped by the island, the electron 
is rapidly energized by the strong (^0.1VaBo/c) reconnec- 
tion electric field. We confirmed that the work done by the 
E z component, i.e. (-e) J v z E z dt, dominates over the work 
done by the E x and E y components, i.e. (-e) J v x E x dt and 
(-e) / v y Ey dt. 

Now, our interest goes to the next period during which 
the energy increases almost constantly until Q c it^ll0 (Fig- 
ure 3a). The energy increment is about 0.3m e c 2 and the 
displacement along the z direction is quite large, Az ~ 
llOdi. What is important here is that the island is expand- 
ing rather than contracting during the energization (Fig- 
ure 2). Therefore, we must explore energization mechanism 
other than the contracting island mechanism [Drake et al., 
2006]. After about Q c it^ll0 when the island reaches the 
edge of the diffusion region, the energy remains almost con- 
stant at e ^1.5m e c 2 . 

The source of energy can again be studied from the pro- 
file of the electric field felt by the particle (Figure 3b). The 
x and y components of the electric field (to be discussed in 
more detail below) are oscillating so that the net energy gain 
from these components should be small. Here, we focus on 
the z component of the electric field, E z , which is almost 
equivalent to the reconnect ion electric field. It remains, on 
average, constant at E z ^0.025VaBo/c until Q c it^ll0 and 
then decreases to zero afterwards. Let us briefly check if this 
E z can explain the energy increment Ae ^0.3m e c 2 . Gener- 
ally, the energy increase Ae is estimated as 



Ae 



Oce \ / cE \ /Az\ 

CJpeJ \VaBq) U J 



(1) 



where Q C e/cJ pe =l/3 in our simulation setup. For the 
particle we analyzed, it felt reconnect ion electric field 
E z ~0.025VaBq/c and was displaced Az ^llOdi during the 
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Figure 2. The trajectory of the most energetic electron during (a) 86. 7< Q c it <90.7, (b) 90. 7< 
n ci t <93.3, (c) 95.0< n ci t <95.6, (d) 103.7< Q ci t <104.3, (e) 112.4< Q ci t <113.0. The background 
images are the electron current density at (a) £l c it=88.7 (b) Q c it=92.0, (c) Q c it=95.3, (c) Q c it=104.0, 
and (d) Q c it=112.7. In (e), the whole trajectory (0< Q c it <115) is shown in a format of y versus e. Each 
red mark shows the time of the background image. 



period 100< Q C it <110 so that the energy increment should 
be As ^0.3m e c 2 in agreement with the above measured 
value. From these arguments, we conclude that the elec- 
tron is energized by the reconnection electric field E z while 
being trapped within the secondary island. 

One may point out a possibility of adiabatic heating be- 
cause the magnetic field magnitude shows gradual increase 
during the energization (Figure 3c). However, this is not 
the case. The spatial scale of the particle motion is compa- 
rable to the island size so that the trajectory is decoupled 
from the magnetic field lines, indicating that the electron 
is highly non-adiabatic (Figure 2c-e). We also confirmed 
that the magnetic moment \i of this particle is not constant 
during the energization (not shown). 

The main feature of the secondary magnetic island within 
a diffusion region is that it contains an electrostatic potential 
well. In the diffusion region, ions are unmagnetized while 
electrons are more magnetized so that the in-plane, polar- 
ization electric field E p is generated and converge toward 
the center of the magnetic island. 

Figure 4a and Figure 4b respectively shows the poten- 
tial and the electromagnetic field structure at and around 
the secondary magnetic island shown in Figure 2c. It is ev- 
ident that the potential well exists in the magnetic island 
and that its depth is of the order of the kinetic energy of the 
reconnection outflow <I> ^miV^- The deep potential well is 
represented by the large, bipolar signature in the E y profile. 
A small bipolar signature also appears in the otherwise flat 
profile of E z because of the motion of the entire island with 
the speed ^0.3Va directed toward the negative ^/-direction. 
This structure led to the oscillatory profile of the E z felt by 
the particle (Figure 2f). Note again that the net energy gain 
is positive because of the E z offset due to the reconnection 
electric field. 



3. Discussion 

We found an energization mechanism by which electrons 
are trapped by a small scale magnetic island so that they 
are continuously energized by the reconnection electric field 
within the diffusion region. This is similar to the 'surfing' 
mechanism by which electrons are accelerated by the recon- 
nection electric field while 'surfing' along the layer of the 
in-plane, polarization electric field [Hoshino, 2005]. In or- 
der to discriminate the two mechanisms, we term the newly 
found mechanism as the 'island surfing'. Let us now con- 
sider the 'island surfing' in a similar way as discussed by 
Hoshino [2005]. As an example, we consider an electron 
that is located along the current sheet center (x—xl) but 
slightly away from the island center (Figure 4c). From the 
symmetry of the electromagnetic fields at and around the 
magnetic islands, the following arguments can be applied to 
other electrons without losing generality. The ^/-component 
of the equation of motion reads 



dp y , 
dt [ ~ 



■e)E y + 



(-! 



(v z B x - V X Bg 



(2) 



where p y =m e Vy is the electron momentum, E y on the right- 
hand side can be replaced by the polarization electric field 
E p produced within the magnetic island, and the guide 
field B g =0 in the present case. If the Lorentz force Fb—{- 
e/c)v z B x is larger than the electric force FE—{—e)E Pl the 
electron is trapped and will gain energy from the reconnec- 
tion electric field E z . From the trapping condition, we get 
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indicating that the 'island surfing' mechanism requires a pre- 
acceleration to be kept trapped within an island. In general, 
the polarization electric field E p is expressed as 



where / is the thickness of the potential structure and may 
be estimated as an intermediate scale between ion inertia 
and the electron inertia length, / = c/^/cjpeCJpi [Hoshino, 
2005]. In the present case, / is measured from the radius of 
the secondary island and is of the order of the ion inertial 
length. This is not unrealistic because the size (diameter) of 
the magnetic islands observed by the Cluster spacecraft was 
about two times the ion inertia length [Chen et al., 2008]. 
For B x , we can assume to be a minor factor smaller than 
Bo from our simulation result. Then, the minimum value 
of \v z \ is a minor factor larger than Va- This is not a se- 
rious condition given the fact that electrons have a great 
chance of pre-acceleration at the X-line as we have already 
seen in Figure 2a. Note also that, in the region away from 
the X-line but within the diff usion re gion, electrons can be 
bulk accelerated up to VAe= \J tth / m e VA [e.g. Hoshino et al., 
2001]. 

An accelerated electron would be kept trapped as long as 
its gyro-radius tl = p e /eBi remains smaller than the size 
of the structure R where p e is the electron momentum and 
Bi is the local magnetic field magnitude felt by the particle. 
Since magnetic field lines of the inflowing plasma constantly 
accumulate on the island, we may assume that the size in- 
creases with the rate a r VA where a r is the instant recon- 
nect ion rate. For the initial size of the magnetic island Ro 
and the time from the start of the island expansion At, R 
= Ro + avVAAt. Then, the trapping condition tl < R can 
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be rewritten as 

1±. < (9-) (*» + aAi At) (5) 

m e c y m e \B J \uj pe J \ di J 

but Ro > di, cjpe > ^ci and Bi/B ~1 in our simulation and 
in the magnetotail as well so that the maximum velocity 
would be 

^z,max ~ C (6) 

in the non-relativistic regime, indicating 'unlimited' electron 
acceleration by the 'island surfing' mechanism (If the elec- 
tron is actually energized up to a relativistic energy, equa- 
tion (5) gives the upper limit). This efficient acceleration 
is an advantage of the closed field line geometry of the sec- 
ondary islands. 

Note, however, that the 'island surfing' works only for the 
island located within the diffusion region entirely covered by 
the reconnection electric field (Figure 4c). Any other island 
outside the diffusion region lacks the reconnection electric 
field and has to generate motional electric field to acceler- 
ate electrons. This can only be achieved by a contracting 
motion (Figure 4d). A bulk motion of the entire island does 
not accelerate electrons because it generates a positive out- 
of-plane electric field at one end and a negative out-of-plane 
electric field at the other end, resulting in a small, net energy 
gain. 

In our simulation, the electron in Figure 2 and 3 
has already reached e/m e c 2 =l with the displacement of 
Az/di ^150 during the initial interaction of the newly 
emerging island. Note that we used an unrealistically high 
initial temperature e t h=9.2x 10 -3 m e c 2 ~ 4.7 keV due to 
limited computational resources. If we normalize simulated 
particle energies by e t h which is ^ 1 keV in the magnetotail 
and assume di=500 km, the simulation indicates electron en- 
ergization up to 110 keV with the dawnward displacement 
^12Re- Given the fact that the diameter of the magne- 
totail is ^40Re, the initial interaction with the emerging 
secondary island may explain energetic (<127 keV) elec- 
trons observed by the Cluster spacecraft that passed through 
small scale magnetic islands [Chen et al., 2008]. 

The main part of the 'island surfing' may work in much 
larger scale current sheets such as those found in the solar 
atmosphere. In this respect, we would also like to point 
out that magnetic reconnection in a spatially large current 
sheet (in the ^/-direction) would yield not just one magnetic 
island but a chain of islands [Loureiro et a., 2007] so that a 
multitude of 'island surfing' builds up to create a significant 
amount of energetic electrons. In addition, if an external 
force existed to drive the magnetic reconnection, the recon- 
nection electric field increases [Hoshino, 2005] and the fre- 
quency of secondary island generation also increases [Wan 
and Lapenta , 2008] . Such a driven configuration is favorable 
for more energization by the 'island surfing'. 

We did not consider a guide field, but our preliminary 
simulation with a guide field (B^=0.3Bo) also showed effi- 
cient electron energization within a secondary magnetic is- 
land, although its efficiency is still under investigation. 

Finally, we would like to point out that the 'island surf- 
ing' is not just a matter of energization of a few unusual 
electrons. It would give an insight toward the understand- 
ing of the roles of secondary islands in the context of energy 
disspation during magnetic reconnection. 
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Figure 3. Time profiles of (a) the energy and the dis- 
placement Az, (b) the electric field and (c) the magnetic 
field felt by the particle during the trapping motion. In 
order to eliminate the particle noise, the electric field data 
have been smoothed by box average with the box size of 
Q ci At=2.7. 
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Figure 4. (a) The electric potential at Q c it=104. The 
contour shows the magnetic field lines, (b) A ID cut 
through x — xl of B x /B , cE^/VaB , and cE^/VaB . 
(c) A schematic illustrations of a secondary island in the 
diffusion region (DR) and the force balance of an elec- 
tron, (d) An illustration of a contracting island. 



